Magnon Dispersion and Anisotropies in SrCu2 (303)2 
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We study the dispersion of the magnons (triplet states) in SrCu2 (603)2 including all symmetry- 
allowed Dzyaloshinskii-Moriya interactions. We can reduce the complexity of the general Hamilto- 
nian to a simpler form by appropriate rotations of the spin operators. The resulting Hamiltonian is 
studied by both perturbation theory and exact numerical diagonalization on a 32-site cluster. We 
argue that the dispersion is dominated by Dzyaloshinskii-Moriya interactions. We point out which 
combinations of these anisotropies affect the dispersion to linear-order, and extract their magnitudes. 



I. INTRODUCTION 

Strontium copper boron oxide (SrCu2(B03)2) is a two- 
dimensional antiferromagnet with no long-range mag- 
netic order. ^'^ Heisenberg interactions between the S — 
1/2 magnetic moments are strong and the geometry is 
such that the planes of spin dimers are coupled in a way 
equivalent to those of the Shastry-Sutherland model.-il 
This model had been proposed on purely theoretical 
grounds, and is remarkable in that while strongly inter- 
acting, it has an exact dimer ground state, fully isotropic 
in spin-space, up to a critical value of the coupling. - It is 
not integrable, however, and neither the excited states, 
nor their dispersions are known exactly. This motivated 
many studies of the dynamics of this model. The ex- 
istence of a compound which is well modeled by the 
Shastry-Sutherland Hamiltonian, and what is more, is in 
an intermediate range of coupling, not far from a quan- 
tum critical point, is then of great interest. The spin ex- 
citations of SrCu2(B03)2 have been studied extensivel 
by a variety of experimental techniques (see e.g. 0| 
and references below). An immediate question is how to 
extract the couplings and to give a precise answer to the 
question: to what extent SrCu2 (603)2 can be described 
by the Shastry-Sutherland model? 

The first observation is that while the Shastry- 
Sutherland model and its ground state are rotation- 
ally invariant, SrCu2 (603)2 shows departures from spin 
isotropy. Electron spin resonance (ESR),^ far- infrared 
spectroscopy of forbidden transitions,* neutron inelastic 
scattering)^ or latter Raman scatteringi^ have shown the 
splitting of the spin triplet excitations at q = and 
anisotropic behavior with respect to the direction of an 
external magnetic field. This was explained^ as originat- 
ing from a Dzyaloshinskii-Moriya interaction (a correc- 
tion of spin-orbit origin )ii whose characteristic vector is 
oriented perpendicular to the copper planes, along the c- 
direction of the crystal. While spin-orbit interactions al- 
ways introduce anisotropies at some energy scale, it was 
a surprise that in a frustrated model they could dom- 
inate the dispersion which generally arises from larger, 
isotropic, couplings. More recently, other components of 
the vectors were argued to be required to explain further 



experimental findings: another splitting was observed at 
q = (tt, 0)f^ii3 and was subsequently associated with 
in-plane components of the Dzyaloshinskii-Moriya inter- 
actions. In addition, there is an avoided level crossing 
when the first triplet state is about to cross the ground 
state for magnetic fields applied along c. This is not com- 
patible with a single-axis anisotropy along the same di- 
rection, and requires additional forms of anisotropy. The 
level anti-crossing was argued to be due to the nearest 
neighbor Dzyaloshinskii-Moriya interactiou j^^d^d^ which 
was considered as a possible explanation for the X-band 
ESR as welljiii All these components are indeed allowed 
by the symmetry of the crystal structure in the low 
temperature phase of SrCu2 (603)2?^'^'^ whose copper 
planes are slightly buckled. Were it not buckled (as in 
the high temperature phasei^), these components would 
have vanished. It was for this reason that they were 
neglected in the original interpretation of the neutron in- 
elastic scattering!^ Nonetheless, because the symmetry 
is slightly broken these interactions are expected to be 
present and define a smaller energy scale. 

A precise model for SrCu2 (603)2 must therefore in- 
clude, in addition to the larger Heisenberg couplings, a 
complex pattern of Dzyaloshinskii-Moriya interactions. 
Since their relative strengths are not precisely known, it 
is difficult to establish a hierarchy among them, other 
than the probable dominance of the c axis component. 

In this paper, we argue, on the basis of exact trans- 
formations and numerical diagonalizations, that we can 
simplify the model and determine the parameters of 
the transformed model by studying the dispersion of 
the triplet states. The dispersion of the triplet states 
of the Shastry-Sutherland model, taken alone, is very 
small, and this is why smaller interactions are rele- 
vant, and even turn out to be dominant?^ It is natu- 
ral to consider Dzyaloshinskii-Moriya interactions since 
they are linear in the ratio of the spin-orbit coupling 
to the crystal-field splitting. A, which is estimated from 
the (^-factor, i.e. to be about 0.1. They are, 

however, peculiar, in that they are usually expected to 
have little effect on the spectrum, at most A^.^^ It has 
even been claimed that Dzyaloshinskii-Moriya interac- 
tions and second-order symmetric anisotropic couplings 
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conspire to restore the rotational invariance of the spin 
excitation spectrum. While this can be true for single- 
band models in certain geometries, here the frustration of 
the lattice makes such arguments inapplicable, as we shall 
see. We can in fact classify the Dzyaloshinskii-Moriya 
components into reducible components (that can be re- 
duced to effects) and irreducible components that have 
effects on the dispersion linear in the spin-orbit coupling. 
The latter thus defines a larger energy scale which will 
be the focus of our discussion. 

The problem then is to calculate the spin excitation 
spectrum reliably. If only isotropic interactions are con- 
sidered, SrCu2 (603)2 is very close in parameter space 
(J'/ J — 0.62) to a quantum critical point, ^ {J' /J — 
0.68), which separates the dimer phase (which is known 
exactly) to a phase that is not known exactly but may 
be quadrumerized.— i^i^ Perturbative approaches to the 
excitations may therefore turn out to be inaccurate. For 
instance, the splitting of the q = triplet energy is given 
hy S — D'j^g{J' / J), which is linear in and involves 
a function g{J'/J) with g{0) = 4, but g{J'/J) = 2.0 
for J'/J = 0.63^. This is difficult to calculate by per- 
turbative techniques because of the proximity with the 
quantum critical point. Therefore, in order to explain 
the dispersion of the first triplet states in SrCu2 (603)2, 
we need not only to include the relevant Dzyaloshinskii- 
Moriya interactions, but also to go beyond the perturba- 
tive techniques used earlier for the dispersioni^ i^^i^^ In- 
deed, even carried at high order in perturbation theory, 
such calculations may tend to overestimate the disper- 
sion, as we shall see. For these reasons, we have carried 
out an exact numerical study of the model, including the 
Dzyaloshinskii-Moriya interactions. 

In addition to symmetry-allowed components of the 
Dzyaloshinskii-Moriya vectors, other components, that 
are not allowed by the measured crystal symmetry were 
invoked2ii2^ as possible explanations for the specific heat 
and the ESR "forbidden" transition. We shall not con- 
sider them here because, on one hand, no distortions of 
the crystal structure that would allow them have been re- 
ported so farjiSi^ and electric-dipole-transitions provide 
a symmetry-preserving alternative to explain forbidden 
transitions. ^'^^ Although it is not clear whether ESR are 
of magnetic-dipole or electric-dipole nature, it has been 
shown that the forbidden transitions observed by far in- 
frared spectroscopy were clearly electric-dipolcj^S Even 
if the ESR transitions were magnetic-dipole, they could 
possibly be explained by finite-temperature effectsi^ 
These components, not allowed by the symmetry, seem 
therefore unnecessary at our current state of knowledge. 

Those components which are allowed by symme- 
try, however, are most probably important to under- 
stand refined experiments, such as the high-field plateau 
phases^^i^ where typical condensation energies are much 
smaller than the isotropic couplings. For instance, the 
1/8 plateau width is of order of 1 Kelvin and it is ques- 
tionable whether Heisenberg models alone capture the 
essential superstructure of that phase. This provides fur- 



ther motivation to quantify to what extent the isotropic 
model applies to SrCu2 (603)2, and in what circum- 
stances the anisotropic couplings can be safely neglected. 

The paper is organized as follows. In section [Til we 
give the dispersion of the first-triplet excitations for the 
Heisenberg Shastry-Sutherland model. In section IIII Al 
we present the anisotropic model for SrCu2 (603)2 with 
Dzyaloshinskii-Moriya interactions, and in piI6p map 
it onto a simpler model. We then calculate the magnon 
dispersion both perturbatively (section llV)) . and by exact 
diagonalization of finite lattices (section |V| . We discuss 
experimental results in section IVII and extract the cou- 
plings that previously could not be quantified by com- 
paring exact spectra with experimental results. 

II. DISPERSION OF THE LOWEST TRIPLET 
EXCITATION OF THE SHASTRY-SUTHERLAND 
MODEL 

We first consider the Shastry-Sutherland modelS. (with 
no Dzyaloshinskii-Moriya interaction), as originally used^ 
to describe SrCu2 (603)2: 

i7 = ^ JS,.Sj +^ J'S^.Sj (1) 

nn nnn 

where nn stands for nearest neighbors (intra-dimer) and 
nnn for next nearest neighbors (inter-dimer) . The lattice 
of couplings is as shown in Figure [U nearest neighbour 
couplings are on the diagonals in every second square 
of the lattice of next-nearest-neighbour couplings (in the 
real lattice the angles are such that the diagonals are 
shorter). The strengths of the Heisenberg couplings 
are J = 85 Kelvin for the nearest-neighbor coupling, 
and J' = 54 for the next-nearest-neighbor coupling, as 
determined^ from susceptibility measurements. An in- 
terplane coupling is present, but it is small (J" = 8 K) 
and frustrated.'' We therefore restrict ourselves to a two- 
dimensional version of the model and neglect J". The 
ratio J'/J ^ 0.62 — 0.63 was confirmed to be consistent 
with the ratio of the first Raman-active singlet energy 
to the first triplet energy^S. a measure independent of the 
magnetic susceptibility. 

We calculate a few low-lying states of the Shastry- 
Sutherland model [Hamiltonian ^] on a 32-site cluster 
using exact diagonalization. This cluster is particularly 
attractive compared to smaller ones not only because 
it has smaller finite size effects, but also because it al- 
lows access to the important symmetry points (7r/2, 7r/2), 
(tt, tt) and (7r,0) [34ii]- The results are given in Fig. [T] for 
J'/J = 0.62. 

The gap between the ground state and the triplet state 
is 0.442 J for the 32-site cluster and essentially indepen- 
dent of the system-size as shown in table ID Yet the ab- 
sence of finite-size effects on the gap does not mean that 
all energies have reached their thermodynamic limit. In 
fact the dispersion of the 16-site cluster is much flatter 
than that of the 20 or 32-site clusters. It has appeared 
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FIG. 1: (color online). Dispersion relation of the first triplet 
excitation of the Heisenberg Shastry-Sutherland model (dia- 
monds), as calculated by exact numerical diagonalization on 
a 32-site cluster, for J'/ J = 0.62. The vertical axis shows the 
excitation energy of those states with respect to the ground 
state. The degeneracy of the two triplets at (7r/2,7r/2) is 
slightly lifted, reflecting the two dimers per unit cell. Singlet 
states nearby are shown as open circles. 



before (but for the simple square lattice) that the 16- 
site cluster can be quite special compared with larger 
clusters.-^ Comparison of the bandwidths of the 20-site 
and 32-site clusters is also difficult because the recip- 
rocal points are not identical. In particular the band- 
width of the 20-site cluster is inherently smaller because 
of the absence of the (tt, tt) reciprocal point: the only non- 
zero point in reciprocal space is at (^, ^). Therefore 
while finite-size effects are definitely smaller than the ta- 
ble H] might suggest it is not possible to quantify them at 
present. In absolute values, the bandwidth of the 32-site 
cluster remains very small, ^ 0.023 J, i.e. about 2 Kelvin. 
We note that it is smaller than that of the perturbative 
calculation,—!^ although this may be explained in part 
by finite-size effects, but is also qualitatively different in 
shape. Consistent with the fact that there are two dimers 
per unit-cell, we have in fact not one but two triplet states 
which are degenerate at high symmetry points of recip- 
rocal space. The degeneracies are slightly lifted for inter- 
mediate q-values. This splitting was shown to occur at 
higher order, (J'/J)^° in perturbation theoryi^ The gap 
and especially the dispersion relation are in fact much 
more affected by the smaller anisotropic interactions, as 
will be shown in the next section. 



III. MODEL AND EXACT MAPPING 

A. The general anisotropic model 

We now consider a model with Dzyaloshinskii-Moriya 
interactions, based on the symmetries of the low tem- 



cluster size 


spin gap 


bandwidth 


16 


0.460 125 J 


0.003 15J 


20 


0.451 441J 


0.015 98J* 


32 


0.442 026 J 


0.023 16J 



TABLE I: The spin gap and bandwidth of the first triplet 
excitation in the Shastry-Sutherland model for different clus- 
ter sizes. *The bandwidth of the 20-site cluster cannot be 
directly compared with the others because of the absence of 
the {n,n) reciprocal point. 

perature phase li^ii^i^S The model contains isotropic 
couplings,* and both inter-dimer— and intra-dimer 
Dzyaloshinskii-Moriya interactions It reads: 

H - ;^[./S,.S, +D,,.(S, xSj)] 

nn 

+ 5][j'S..S,+D^.(S,xS,)] (2) 

nnn 

where nn stands for nearest neighbors (intra-dimer) and 
nnn for next-nearest neighbors (inter-dimer). The dif- 
ferent components of the vectors Dij and D^^ are given 
in Fig. [5) They are obtained by using transformation 




FIG. 2; (color online). Shastry-Sutherland lattice with 
Dzyaloshinskii-Moriya interactions allowed by the symmetries 
of SrCu2 (60^)2. Colored (short gray) arrows are the com- 
ponents of the Dzyaloshinskii-Moriya vectors for each bond 
(long arrows show the prescription for the order of the op- 
erators, a black arrow from i to j indicates that we should 
take D • (Si x Sj) with the given D). mi and 7712 are mirror 
planes, but the (ab) plane is not a mirror plane in the low- 
temperature phase. This allows for all in-plane components 
of D. 

properties of pseudo-vectors and the symmetries of the 
crystal structure For instance, when a bond connecting 
two spins belongs to a mirror plane, the Dzyaloshinskii- 
Moriya vector must be perpendicular to it. In the high 
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temperature phase (T > 400 K)^ the {ab) plane is a 
mirror plane and the allowed components (only the inter- 
dimer ones) must therefore be parallel to the crystallo- 
graphic c-axis. This is the D^^ — D'^Gz (e^ is the unit 
vector along z — c) we considered before. In the low 
temperature phase, the {ab) plane is slightly buckled, 
and some other in-plane components are allowed, but are 
expected to be smaller The remaining toi, m2 mir- 
ror planes allow us to determine the complete pattern 
of Dzyaloshinskii-Moriya vectors. In-plane components 
of Ti'^j can be decomposed into staggered Dy ^ and non- 
staggered _D|j orthogonal components. We obtain 

D = D\ ^,e,,,, + D\ „,e„,,,, + D'^b, (3) 

where the unit vectors are, respectively, es,ij = ±ej, and 
Bns^ij = ±ex (for dimers A) or e^^ij = zbe^^ and ens,ij = 
zbCy (dimers B) (Fig. [2]). The intra-dimer vector has 
only in-plane components, 

D,, = De,, (4) 

with the unit-vector e^j = Gj, (dimers A) or Gij — Qx 
(dimers B). Note that because the interaction is a vec- 
tor product of spins, we have to be careful with the 
order of the operators in the definitions. The order 
is defined in Fig. [5] for each bond by the black ar- 
rows. For instance, for the bond 1 — > 2, we define 
'D'i2 — s^y + -^11 ns^x ~ D'j^e^- The strengths of three 
of these four couplings, D, Dy ^, Dy are unknown. In 
the next paragraph, we show how to simplify the Hamil- 
tonian Q such as to eliminate two of these couplings. 
Eventually only one parameter remains to be determined 
and we will determine it with the help of the magnon dis- 
persion at finite wave-vector. 



B. Mapping 

We map the Hamiltonian ^ onto a simpler Hamilto- 
nian with no intra-dimer and no uniform Dzyaloshinskii- 
Moriya interactions, by means of rotations^^ of spin op- 
erators. Contrary to previous situations where this ap- 
proach has been appliedj^ii^^ it turns out here that we 
cannot eliminate all components, because of the frustra- 
tion of the lattice but we can nonetheless simplify the 
model. For this, we perform two sorts of rotations. First, 
we rotate the two spin operators of the dimers in oppo- 
site directions, and apply the same operation to all unit 
cells. Second, we rotate the spin operators from dimcr 
to dimer by applying the same operation to the two spin 
operators of the same dimer. The first transformation 
uses local rotations (the same for all unit-cells) for the 
two spins i and j of the dimer i, j, that are defined by: 

S', =7^(0/2,ey■)S, (5) 

9 9 6 

= (1 - cos -)[ey.Si]ey -I- cos -S^ - sin -S-^ x Gy (6) 



with tan 9 = D/J^ where D and g^ are defined by ([4]). 
The second spin j of the same dimer is rotated in the 
opposite sense {9 —S). For the second dimer of the 
unit-cell, the spins are rotated similarly, about an axis 
that is perpendicular to that, according to the local Dij 
(see Fig.EJ. 

We now restrict consideration to terms that are lin- 
ear in the Dzyaloshinskii-Moriya coupling strength in the 
transformation ([6]), because we are looking for linear ef- 
fects in the dispersion of excitations. In fact the trans- 
formation ^ does produce some anisotropic exchange of 
order / which we neglect because (i) we have not 
taken into account real anisotropic couplings which occur 
at the same order of magnitude, and more importantly, 
(ii) the spectrum is fully dominated by effects that are 
linear in D, as we shall show. 

In terms of the primed spins, the Hamiltonian ([2]) is 
mapped onto (up to terms of order / J): 

H^Y. ^S:-S; + E [J'^'^■^'3 + D:,.(S^ X S;)] (7) 

nn nnn 

where the D^^- are now modified: the staggered c- 
component is unchanged, whereas the strengths of the 
in-plane components (either staggered or uniform) are 
modified according to 

Dins ^ Dlns + ^jD (8) 

D'i_ ^ D'^ (10) 

This is an exact transformation that is valid at all order 
in J' / J. As a second step, we now proceed to eliminate 
the non-staggered component of D-^ , which is defined by 
its strength Dy (eq. ([8])) and a vector direction either 
along X 01 y (see Fig. [2]). This is made possible because 
the sum up to zero along any closed loop of the 

lattice. When going around closed loops, the rotation 
angles are fixed by the Dzyaloshinskii-Moriya strength 
and sign; therefore there is a compatibility condition for 
the last spiniS Consider for instance the closed loop (2 — 
8 — 5 — 3 — 2) in Fig. [21 the sum of the Dy is zero along 
this loop. Note, however, that the sum of the D'j^ does 
not vanish, nor does the sum of the Dy ^ along (5 — 7 — 
6 — 5), for instance, so these cannot be eliminated. We 
now rotate the operators accordingly, and restrict to the 
first-order terms in £)|| 

S, = ni9,,ey)n{iP,,ex)S[ (11) 

In this expression, 9i = (I?|| J')Ri • e^, and ipi — 
i-D'w ns/"^')-^* ■ ^i" where e^^y are unit vectors along the 
X and y directions. is the position of the dimer to 
which the spin i belongs, in units of the inter-dimer spac- 
ing. The two spins of the same dimer are rotated in 
the same way, so as not to generate other intra-dimer 
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Dzyaloshinskii-Moriya interactions. The final Hamilto- 
nian becomes: 



T>,,,=Dne 



D\\ = D'n 



2J 



(12) 

(13) 
(14) 



where e^^ij is the unit-vector along the inter-dimer stag- 
gered Dzyaloshinskii-Moriya vector, defined in Fig. [2l 
The transformed model (fT2|) is simpler and can be used 
instead of the original Hamiltonian 12]), as far as ener- 
gies are concerned. The remaining Dzyaloshinskii-Moriya 
components (1131) which we call irreducible cannot be 
eliminated by rotations. Geometrically this is because 
they do not sum to zero when turning around closed loops 
of the Shastry-Sutherland lattice. As a consequence, lin- 
ear effects are present in the spectrum. 

We shall proceed with the Hamiltonian (I12|) since the 
two Hamiltonians have exactly the same energy spectrum 
up to terms of order -D^/J, D'^^^^j J . Note that this is 
valid only at zero magnetic-field: at finite fields along 
0, for instance, the first transformation creates a stag- 
gered transverse field^^ (which must already be present 
according to the local environment surrounding the cop- 
per ionsi^) and the second transformation creates rotat- 
ing fields. 



IV. PERTURBATIVE CALCULATION OF THE 
EXCITATION SPECTRUM 



We calculate the magnon dispersion and dynamical 
structure factor of model p^ . by first-order perturba- 
tion theory in the Dzyaloshinskii-Moriya couplings and 
J'/Jj£vl2, Since the latter is not small for SrCu2 (603)2 
(J'/ J = 0.62 — 0.63), we will supplement this approxima- 
tion with exact numerical results in section |V] but per- 
turbation theory already gives some qualitative insights 
into the problem. 

We start with the exact dimer ground state, which is 
a product of singlets on each dimer. We consider the fol- 
lowing one-particle excited states, which consist of pro- 
moting a dimer into a triplet state with 5^ — 0, ±1 quan- 
tum number. They are the eigenstates for zero J' and _D, 
and we use first-order perturbation theory in J' j J , D\\ 
and D_L . Because there are two dimers per unit-cell, there 
are two such triplet states depending on which dimer A 
or B is in the triplet state. When is added, the total 
rotation invariance is broken but 5^ is still a good quan- 
tum number thanks to the rotation about c. In this case, 
the eigenstates within this subspace of states are given 

\s\ ±> = ^ q> A) ± i\s\ q, b)) (15) 

with S'- = 0, ±1 the z component of the total spin of 
a dimer, and |S'^,q, A) is the Fourier transform of the 



triplet state on the dimer A of unit-cell z, ^i^A). Now 
in addition we consider the effect of inter-dimer stag- 
gered components of the Dzyaloshinskii-Moriya vectors 
i.e. model (jl2p . Dy breaks the rotation invariance about 
c and therefore mixes various S"^ components of (jlSp . 
The Hamiltonian within the subspace of splits into 
two blocks. The first block for (| + 1, q, |0, q, -), | - 
1 , q, +) ) reads 



J + ^^/(q) -i^||5(q) 
z^yg(q)* J 

-^^||.g(q)* 





*^ii5(q) 

J-^±/(q) 



(16) 



where we have /(q) = 2 cos((7a/2) cos((7b/2) and ^(q) — 
- (sin[(gfa + gh)/2] -I- isin[((7a - 9b)/2]) /\/2, q= {qa,qb)- 
The lattice spacing is taken such that a — 1. The second 
block for the three other states, is identical up to the 
change (l)|| , D±) , D±). The six resulting triplet 

states are denoted by and have energies (twice de- 
generate because of identical blocks) : 



i?±(q) = J±y'^i/(q)2 + 2^||5(q)P^ J±f 



E"{q) = J 



(17) 



The dispersion relations of "triplet" states pT)) are given 
in Fig. [3l They are valid only at first-order in J'/ J and 
higher-order corrections are not included here. We shall 
see below how the overall shape and degeneracies are 
modified when one goes beyond perturbation theory. For 
the moment, this defines two splittings at q = (0,0) and 
q = (tt, 0) that are given by: 



5(0,0) 



= 4Dj_ = AD'^ 

= 2\/2i)|| = 2\/2 [d\ 



2J' 



D 



(18) 
(19) 



where the right-hand sides give the expressions of the 
splittings in terms of the original couphngs of Fig. ^ 
according to the transformations © and PH]) . It is in- 
teresting to note that measuring these splittings at dif- 
ferent wave- vectors allows, in principle, separation of the 
two couplings. This is because of the different staggering 
patterns of the Dzyaloshinskii-Moriya vectors in space, 
which is itself a consequence of the crystal symmetry. 

In order to compare with neutron inelastic scattering 
at finite wave- vectors, we have computed the neutron 
cross-section, given by^ 



\dLjdQ J 

a/3 ^ ^ 



5"(q)<5(c.-i?"(q)) (20) 



(0|5^|t^')(i^ 15^10) (21) 



Note that 5q is, strictly speaking, the real spin opera- 
tor, not that obtained after the transformations are per- 
formed, S'q. We can in fact replace S'q by S'^ because it 
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FIG. 3: {J' / J 0) Dispersion relation and dynamical struc- 
ture factor (grey scale) of the first triplet excitations in the 
limit of small J' /J. Note the overall dimer form factor Mq is 
not included in the picture. For (7r,0), the middle mode has 
no intensity. For the clarity of the picture, D_i_ = O.IJ, and 
D\\ = 0.05J were used. 



implies small intensity corrections, of order D/J, which 
we shall discuss later in the paper. The cross-section 
involves calculating all matrix elements. Some do not 
contribute for special directions of the reciprocal vector, 
such as for those along the reciprocal line q ~ {q, 0) 
(q parallel to a gives geometrical factors Qa/q^ = 1, 



= 0), and only (t; 
tribute. For q = (g, q) and 
have the form 



q,5^|0) and (i™|5^|0) con 



5±(q) = 



1 , / 2iJ||l.9(q)l 

2 I 5„ 



((7,0), the intensitiei 
(22 

Mq + • • • (23 



where Mq — sit? q • (5 -I- sin^ q • (5' is the form factor of th( 
unit-cell, 25 and 25' being the vectors along the dimeri 
A and B (Fig.[2|). The dots represent corrections of orde: 
©(Z?,!?!! whereas the intensities are 0(1). We wil 
come back to this when discussing the effect of D an( 
D|| For q = (tt, g), the q dependence comes only fron 
the geometrical factor: 



that the intensity of the middle mode at (tt, g = 0) is zero. 
This is consistent with the presence of only two modes in 
neutron scattering experiments at that wave- vector , 
and allows us to identify these modes as the lower and 
upper modes (as will be confirmed by their energies later 
as well). 



V. EXACT NUMERICAL DIAGONALIZATION 

We calculate a few low lying states of the Hamiltonian 
([T2| on the 32-site cluster using exact diagonalization. 
Note that the term in (fT2|) breaks spin rotational 
symmetry and the dimension of the Hilbert space now 
becomes ~ 5 x 10*, as compare to ~ 7 x 10'' for the 
32-site cluster of the Shastry-Sutherland model. Results 
are shown in Fig. H) The parameters are J' / J = 0.62, 
b±_ = 0.02J and b\\ = 0.005J. It is to be compared 
with Fig. [T] which corresponds to D±^ = = 0. Since 
S"^ is no longer a good quantum number, degeneracies in 
those triplet states in Fig.[T]are lifted. Nevertheless, they 
can still be grouped together according to whether their 
energies increase, decrease or remain unchanged in the 
presence of a weak magnetic field along the c direction, 
which correspond to = +1, —1, and respectively in 
the limit iPy 0. In Fig. 3] states belonging to the same 
group at different q arc joined together by a line. Note 
that the middle state, which corresponds to the ^ 
and indicated by diamonds, is almost unchanged with 
respect to Fig. [TJ But the others acquire a larger disper- 
sion that in any case remains small compared to the gap. 
Next we study in detail the effects of D± and D\\ on 
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The intensities (eqs. (P^ - d^SI) ) are represented in grey 
scale in Fig.[3](Mq is not included). We see, in particular. 



1^"' + ^ 

2 TT^ + q^ ^ 



FIG. 4: (color online) Same as Fig.[T]but for the Hamiltonian 
(fT2)) with J'/J = 0.62, D± = 0.02J and = 0.005J. Lines 
are guides to the eye and join together states which would 
correspond to = 4-1 (triangles), —1 (squares), and (di- 
amonds) in the limit Dy 0. Note that some of them are 
indistinguishable for the upper mode. 
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FIG. 5: (a) Effect of D± on the q = triplet (diamonds) and 
singlet energies (open circles) when = 0, and (b) effect of 
Z)|l on the q = (tt, 0) states when D± — 0.02 J. Numbers are 
calculated using exact diagonalization on the 32-site cluster. 
Lines are linear fittings to triplet energies. 



the splitting of the first triplet state. Fig. El^a) shows the 
splitting in the q = "triplet" state at various D± and 
= 0. We found that a small 0.005 J) has no no- 

ticeable effect on the splitting of the q = states. There- 
fore Fig.[5]Ja) can be regarded as showing the effect of D± 
on the splitting of those states. Fig.jUJb) shows the split- 
ting of the first triplet state at q = (tt, 0) at D± — 0.02 J 
and various -Dy . It is clear that the effect of is to mod- 
ify the dispersion around the reciprocal point (7r,0). At 
(tt, 0), Fig. inib) shows that the energies change linearly 
as function of the coupling, except at very small Dy. In 
this region, a very small gap already exists, and is given 
by D±{J' / J)" (where a > 1). In absolute values we have 
found a gap of 1.6 x lO'^j (for D± = 0.02J). The tiny 
gap is amplified by Z?|| and the effect is perturbative in 



the strength of the Dzyaloshinskii-Moriya interactions. 
Note that the slopes, however, are a function of J' /J 
and are non-perturbative in this coupling. First, for val- 
ues of J'/ J close to the critical point, the separation of 
the components as in and (fTQ]) is no longer exactly 
true, but the effect of -Dy at q = remains small (and 
conversely) in the regime of interest. Second, the prefac- 
tors of the splittings between the upper and lower modes 
are modified compared to the perturbative results. From 
linear fittings as shown in Fig. [51 the splitting at q = 
is found to be ~ 2.08i)j_ instead of 4i)j_ in HH]). This 
compares very favorably with a previous calculation us- 
ing the 20-site cluster^ thus showing that size effects are 
small. Similarly at q = (tt, 0) the splitting is ~ 2.80i)y 
instead of 2\/2i^|| in ([TO)) . Note that the renormalization 
at q = (tt, 0) is in fact very small compared to that at 
q = 0, and solves the puzzle previously noticedfi^ as we 
shall explain below. 

The spectra show other minor features: at q = 0, for 
instance, the lowest state is split into two (Fig. 3]) . This 
is to be expected when there are several anisotropy axes 
and no Kramers degeneracy in a system with an even 
number of spins per unit-cell. We can simply understand 
these additional splittings by perturbation theory: cou- 
plings with higher energy states open gaps which scale as 
Z^jj (which we verified by exact diagonalization). The ex- 
act calculation shows other effects. We cannot make 
quantitative comparison of effects at this order to real 
experiments because anisotropic super-exchange is also 
present at order and are neglected in the present work. 
A^ effects may perhaps be accessible to ESR at q = 0, 
but it might be difficult to disentangle the second-order 
effects of Dzyaloshinskii-Moriya from that of real asym- 
metric couplings. 



VI. EXTRACTING THE COUPLINGS FROM 
EXPERIMENTS 

Triplet states. We now use these results to extract the 
anisotropic couplings from neutron inelastic scattering 
experiments. At q = 0, we have seen that the energy is 
unchanged between the 20-site cluster^ and the present 
32-site cluster, so the estimation of = 0.18 meV is 
unmodified. Regarding the intensity at q = 0, earlier 
perturbative calculation^ compares well with exact nu- 
merical diagonalization on a 20-site cluster ^22, so we have 
not carried out intensity calculations on the 32-site prob- 
lem, but will rely on perturbative results. At q = (tt, 0), 
we have seen indeed that the intensity of the middle mode 
is zero (section ITV)) . This is consistent with the observa- 
tion of vanishing intensity for the middle mode in neutron 
scattering The splitting seen at q = (tt, 0) can there- 
fore be interpreted as the gap between the lower and up- 
per state. As we have seen, the splitting at q = (tt, 0) is 
givenbyi)y/(J7J), with/(0) = 2\/2 and /(0.62) = 2.80 
(which is in fact within 1 % of the perturbative result). 
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This allows to extract Dy ~ 0.07 meV, taking a splitting 
(5(7r, 0) = 0.2 meV— ii^. This is the value we extracted 
before relying on the perturbative expressionji^ We have 
here calculated the renormalization coefficient at (tt, 0), 
and shown that that coefficient is quite different from 
that at (0,0); therefore applying the same renormaliza- 
tion to both splittings^'^ leads to an artificially large in- 
plane Dzyaloshinskii-Moriya coupling. This explains why 
the ratio appeared so large previously. In fact, the cor- 
rect ratio is D\\/D^ ^ 0.4 with D'^ = 0.18 meV. We 
recall that Dy is a linear combination of two interactions 
(see below) , and as such may explain why it is larger than 
expected on the basis of the small buckling of the crystal 
structure. We therefore conclude that from the splittings 
at q = (0,0) and q — (7r,0), we constrain the couplings: 

D'^ = 0.18 meV (26) 
V 

D'w s + 7r-;D = 0.07 meV (27) 

11,5 2J 

These values also explain the dispersion of the excitations 
as shown in Fig. [SI at least for the additional reciprocal 
points available in the 32-site cluster (diamonds). Note 
that the scales of experiment and numerical calculations 
are slightly shifted. The shift could be eliminated by an 
optimized choice of J' / J in the range 0.62-0.63, consis- 
tent with previous estimates. 
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FIG. 6: (color online). Magnon dispersion : exact diagonal- 
ization with Dzyaloshinskii-Moriya interactions (diamonds) 
and experimental data from Ref. (circles). Lines are 
guides to the eyes. Parameters are D±_ — 0.18 meV, and 
= 0.07 meV. J'/ J was kept equal to 0.62. This is re- 
sponsible for the shift in the overall scale of the calculated 
energies. 

How could we, at least in principle, determine indi- 
vidually D, ^ and ^^,7 As we have shown, the 
effect of Dzyaloshinskii-Moriya on the spectrum is large 
only for the two effective components that could not be 
eliminated by our transformations. To determine the 
other components, the effect of an external magnetic field 



may provide more constraints. However, staggered gyro- 
magnetic tensors come into play, that are also unknown. 
Alternatively, we can use neutron scattering intensities. 
First, we calculate the correction to the expressions of 
the intensities (eq. ((22)) or p4)) ) at q = (7r,0), coming 
from the intra-dimer Dzyaloshinskii-Moriya interaction, 
D. Indeed the spin operators have a linear correction 
in D (see eq. ([5])) that we have neglected so far. The 
intensity at q = (vr, 0) (that was zero before) acquires a 
term proportional to (I?/J)^: S'°((7r,0)) = -gjiNq where 
Nq = cos^ q • (5 -I- cos^ q • 5'. So in principle it is possible 
to extract D from measuring the ratio of intensities at 

q = ('^'O)' l°(K0)) = i^tan^ This ratfo is very 

small ~ 0.02(D/J)^, and it therefore seems difficult to 
extract this coupling independently. Second, we consider 
the effect on the intensities of the non-staggered compo- 
nents Dy + J'D/2J, that we had ehminated by the 
second rotation. Since the rotation is performed using 
a constant angle when going from dimer to dimer, the 
spin operators contains operators that are shifted in 
q space, 5'q_|_]^.^, where kp is determined by the coupling 
Z?!! + J D/2J. A consequence is that "ghosts" of the 
original magnon dispersion (as given in Fig. 2]) shifted 
by ±ko must appear in the spectrum. Of course be- 
cause the shift is small, of order A, and the dispersion 
is of the same order, the energies are changed by A^. 
The effect will be difficult to detect by neutron scatter- 
ing using constant wave- vector scans, but could appear in 
constant energy scans. This is reminiscent of magnon dis- 
persions in well ordered systems, such as Ba2CuGe207, 
where the Dzyaloshinskii-Moriya interaction and the he- 
licoidal magnetic structure give rise to three branches in 
the excitation spectrum.^" 

Singlet states. In Figs. |4] and [5l higher and lower sin- 
glet states are also given. Since the rotational symmetry 
is broken, denoting a state as a "singlet" is to be un- 
derstood as its property for vanishing anisotropics. Sin- 
glets are not directly visible by ESR, because the exter- 
nal field is swept at constant wavelengths and only field- 
dependent states are detected. Nonetheless the bend- 
ing of the triplet energy levels tend to signal a level 
anti-crossing with a singlet: experimentally such bend- 
ings were observed and thus singlet energies indirectly 
measured4^ The measured anti-crossings are at 2.67 meV 
(compared to the lowest triplet state at 2.81 mev) and 
3.56 meV (compared to the higher triplet at 3.16 meV).— 
In Fig. [4] the lowest singlet has an energy which is 95% 
of the triplet gap, in good agreement with the experi- 
mental figure. However the second singlet seems to be 
higher in energy compared to that of Fig. 3] We know 
that there are other singlet states at higher energies, and 
that observed by ESR may be one of them. 

Sign of the coupling. We remark that although the 
sign of D'^ is of no consequence on the energy spectrum, 
it interchanges the "handedness" of the wave functions 
of the upper and lower states. It therefore has measur- 
able consequences when polarized light interacts with the 
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tripletsi22. It is interesting to note that the sign of 
extracted by analyzingS2, far-infrared experiments^''^^ is 
opposite to that suggested to explain avoided crossing 
at higher fieldsj^^ii^ The apparent contradiction may be 
resolved because the avoided crossing is not determined 
solely by _D, as suggested earlie r^'^'^^i-^^ . but rather by 
a combination of D'^^ ^ + J'D/2J, D[| „^ + J'D/2J and 
the staggered fields. Quantitative understanding of the 
avoided crossing would need determination of each indi- 
vidual coupling. 

VII. CONCLUSIONS 

We have considered a model for SrCu2(B03)2 with all 
the Dzyaloshinskii-Moriya interactions compatible with 
the crystal structure, i.e. a model with all anisotropic 
couplings linear in the spin-orbit coupling. We have con- 
structed a simpler anisotropic Hamiltonian (|12p with a 
smaller number of couplings, by appropriate mappings, 
which has nonetheless exactly the same spectrum up to 
second-order interactions. The transformation allowed us 
to separate what we term the "reducible" components of 
Dzyaloshinskii-Moriya vectors (the non-staggered com- 
ponents -Djj and part of the intra-dimer interaction D) 
from the "irreducible" (the staggered D^, D'^^ ^ with a 
contribution from the intra-dimer interaction D). By 
definition the latter have effects on the spin excitation 
spectrum that are linear in the strength of the coupling, 
whereas the former have second-order effects. The linear 
effect arises, precisely, because we cannot eliminate these 
components. Had we been able to eliminate all of them, 
we should have taken into account anisotropic symmetric 
exchanges., on an equal footing. In that case, depend- 
ing on the underlying superexchange processes, one may 
recover a fully rotationally invariant spectrum (if single- 
orbital processes are dominant) ^21 or not (when multi- 
orbital processes exist). ^1 In any case, the effects would 
have been at most , and very difficult to extract exper- 
imentally because this would be of order of a few /ieV. In 
contrast, SrCu2 (303)2 represents a real situation where 
the bond frustration leads to an energy scale of order A 
-hundreds of ^eV here- in the spin excitation spectrum. 

We have explained quantitatively the dispersion of 
the lowest triplet states (Fig. [6]) which is dominated by 
Dzyaloshinskii-Moriya interactions (compare Fig. [T] and 
U]). First-order perturbation theory was not sufficient 
to account quantitatively for the dispersion shape (see 
Fig. [3]) but allows us to capture most of the symmetries 
except for the small additional splittings due to the pres- 
ence of two dimers per unit-cell. Measurements of the dis- 
persion under high-pressure^^ should in fact be carefully 
interpreted: the evolution of the bandwidth with pressure 
should primarily reflect the change in the Dzyaloshinskii- 
Moriya couplings and not the change in the ratio J' / J 
that measure the proximity to the critical point. The 



spin gap, where the effect of Dzyaloshinskii-Moriya is 
secondary, may serve as a more sensitive probe for the 
proximity to the critical point. 

The model (jl2p whose parameters have been quanti- 
tatively determined, may serve as a good starting point 
to understand other properties of SrCu2 (603)2 at zero 
field, such as higher energy bound states. These states, as 
observed in particular by ESR,^'^ Raman scattering^^iiS 
far-infrared spectroscopy,^^ and neutron scattering^i^ 
have been analyzed so far within the framework of the 
Shastry-Sutherland uiodel^^^ but should be affected 
by Dzyaloshinskii-Moriya interactions. This may explain 
the splittings observed)^ and resolve the issue of "local- 
ization" versus "delocalization" of these excitations. The 
large dispersion seen first in ref. @ and supported by 
theoretical calculations on the Shastry-Sutherland model 
was revealed latter to be composed of several more local- 
ized excitations.--- This remains to be fully understood 
and model p2)) may help. 

We have focused on zero-field results: finite field effects 
are particularly interesting but involve other unknown 
couplings. If the mappings of section Fill Bl (transforma- 
tions ^ and are made with an external magnetic- 
field, both an effective staggered field)^ and rotating 
fields appear. As a real staggered field is already present 
according to the local symmetry of the copper ions, the 
net effective fields cannot be known precisely. The level 
anti-crossing that was found when the spin gap is about 
to close for a magnetic-field along the c direction, in par- 
ticular, cannot be solely due to the intra-dimer inter- 
action, as previously claimedJ^ It has to be seen as a 
consequence of all these effects: intra and intcr-dimcr 
components, and staggered fields. 

Interesting developments for doped samples of 
SrCu2 (603)2 are underway.^^i^ The doped compounds, 
by breaking translation invariance, may also help to dis- 
entangle the different components of the Dzyaloshinskii- 
Moriya vectors, by studying the induced magnetization 
of the magnetic ions in the vicinity of the dopant, by 
NMR for instance. 
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